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Recent developments in the computation of two-loop master integrals for massive Bhabha scattering are briefly 
reviewed. We apply a method based on expansions of exact Mellin-Barnes representations and evaluate all planar 
four-point master integrals in the approximation of small electron mass at fixed scattering angle for the one-flavor 
case. The same technique is employed to derive and evaluate also all two-loop masters generated by additional 
fermion flavors. The approximation is sufficient for the determination of QED two-loop corrections for Bhabha 
scattering in the kinematics planned to be used for the luminosity determination at the ILC. 



1. Introduction 

Bhabha scattering is the process employed to 
measure the luminosity at electron-positron col- 
liders, because of its clear experimental signature. 
At machines operating at a 1 — 10 GeV centre- 
of-mass energy t/s, the relevant kinematic region 
is that of large-angle Bhabha scattering. Small- 
angle Bhabha scattering, instead, is an invaluable 
luminosity monitor at high-energy colliders in the 
TeV region. 

In order to minimize the luminosity error, a 
precise theoretical computation of radiative cor- 
rections to the Bhabha-scattering cross section 
is required. The electroweak next-to-leading or- 
der (NLO) corrections to Bhabha scattering were 
computed a long time ago in pQ. In recent 
years, studies have been focusing on pure quan- 
tum electrodynamics (QED) contributions be- 
yond the one-loop level. The two-loop virtual 
corrections for massless electrons were obtained 
in [SJ. However, this result was not immediately 
useful since the available Monte Carlo programs 
employ a non- vanishing electron mass m. 



The virtual and real second-order contribu- 
tions to Bhabha scattering, enhanced by factors 
of ln(s/m 2 ) and ln 2 (s/m 2 ) were completed in 
[3141516) . This result was recently improved in 
7 8 9, where the photonic non- logarithmic term 
was evaluated at leading order in the ratio m 2 /s. 
The diagrams with fermion loops remained un- 
covered in this approach. 

An important breakthrough in the field was the 
use of the Laporta-Remiddi algorithm ( [10111] ). 
in order to reduce the Bhabha-scattering cross 
section to a few Master Integrals (Mis). The 
technique of differential equations proved use- 
ful in evaluating several Mis (see i.e. [121131 
114)). The results were represented in terms of 
Harmonic Polylogarithms (HPLs) introduced in 
[15] or of Generalized Harmonic Polylogarithms 
(GPLs) (details in the context of Bhabha scat- 
tering can be found in [16] and in references 
therein). These results led eventually to the ex- 
act result of [17118) for the virtual and real next- 
to-next-to-leading order (NNLO) corrections to 
the Bhabha-scattering cross section involving one 
electron loop. Non-approximated expressions for 



all NNLO contributions, except for double box 
diagrams, can be found in |19|20j . The Mis for 
the loop-by-loop contributions were studied e.g. 
in [H]. 

The complete set of the needed master inte- 
grals is known from [13] . Table [1] reproduces all 
two- loop box master integrals for Nf = 1. Nota- 
tions are exactly those of [13]. The B714m3d2 is, 
e.g., a box MI with 7 internal lines (71), four of 
them being massive (4m), with a higher power of 
one of the numerators (a line being dotted (d)); 
it is one of several such topologies and of several 
ones with dots, so '3d2'. In order to improve the 
Bhabha-scattering theoretical prediction, we in- 
vestigate two classes of NNLO QED corrections. 
In Section [2j we briefly discuss a method based 
on expansion of Mellin-Barnes (MB) representa- 
tions ( |26|27|28I29"] ) and review the results of [23], 
where all planar two-loop box Mis were obtained. 
The non-planar Mis are indicated in Table [1] In 
Section [3] we apply the same method to evaluate 
the Mis arising from diagrams containing heavy 
fcrmions, like muons and tau-leptons; in the fol- 
lowing, we will call them the Nf > 1 contribu- 
tions. Their topologies are shown in Figured] 

The MB-representations are valid for arbitrary 
kinematics. Although their actual evaluations are 
restricted to the high-energy limit (small lepton 
masses at fixed scattering angle), they are well 
suited for practical applications. When dealing 
with Nf > 1 Mis, a second fermion mass M 
is involved. Since our purpose is to evaluate 
the complete QED Bhabha-scattering cross sec- 
tion at high energies, we assume a hierarchy of 
all three scales, namely m -C M -C s, t, where 
t is the usual Mandelstam invariant related to 
the scattering angle. With the summation tech- 
niques described in Section [4] divergent parts 
have been evaluated exactly. Note that the treat- 
ment of hadronic contributions is a separate prob- 
lem, which is better solved by using dispersion 
relations (see [3D]). 




V412M2m B512M2m 

Figure 1 . The topologies of the eight master in- 
tegrals for the heavy-fermion corrections. Bold 
lines represent heavy fermions. 




Figure 2. The diagrams for Bhabha scattering 
with two fermion flavors. Internal fermionic loops 
represent heavy leptons and the other fermion 
lines arc electrons. 



Table 1 

The Nf = 1 four-point master integrals entering the six basic two-loop box diagrams. NP denotes 
non-planar topologies, and references with a dagger give divergent parts only. 
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2. The planar two-loop boxes for Nf = 1 

There are 24 planar two-loop box Mis to be 
determined for the Nf = 1 case, see Table [T] So 
far, we could not determine all of them analyti- 
cally with the exact mass dependence; the status 
is reviewed in |31| . For this reason, we decided 
to treat these Mis uniquely in the approximation 
of small m at fixed scattering angle. The results 
have been published in the mean time [23], so 
we will make here only few introductory remarks 
and show one example. In the list of Mis, we pre- 
ferred to include Feynman integrals without any 
numerators. The pragmatic reason was the inde- 
pendence of their defintion on the internal flow 
of momenta. It is sufficient to indicate the lines 
with dots. Performing explicit calculations, one 
is, of course, faced with the observation that the 
singularities of dotted Mis and those with numer- 
ators are quite different. Evaluating the small 
mass expansions, we preferred in some cases to 



treat instead of dotted integrals those with nu- 
merators. Due to unique algebraic relations be- 
tween all the integrals, there is no principal differ- 
ence in these two approaches, and further details 
are discussed in [23]. Another observation con- 
cerns the MB-representations. In principle, one 
may use the representations for the basic 7-line 
boxes as given e.g. in jJ2J and shrink lines. How- 
ever, when calculating Mis with numerators, ad- 
ditional representations have to be determined. 
We observed further, that it is sometimes not 
evident how to get an effective represenation for 
dotted Mis from more general ones. For the MI 
B512m2 (a 5-liner) we got, by shrinking of lines 
in B714ml (first planar double box)and after ex- 
panding in e, 11 MB-integrals with at most 4 in- 
tegrations. From our direct derivation, we got 4 
integrals, at most 3-dimensional. For the related 
dotted MI B512m2d2, we got from line shrinking 
102 integrals, and by direct derivation only one, 
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Figure 3. The 5-line topology B512m3. The mo- 
mentum distribution has been chosen to make the 
derivation of the MB representation easier. 



again 3-dimcnsional. 

As an example, we reproduce here for the MI 
B512m3(fc2P3), shown in Figure [31 the basic d- 
dimensional MB-represenation: 
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The small mass expansion of the result is: 
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The expression is more complicated than those 
for the planar 7-line Mis concerning both the 
functions appearing as well as the dependence 
on all three Mandelstam variables s,t,u; the lat- 



ter is typical for non-planar diagrams, where the 
B512m3 topology contributes. 

3. The master integrals for the Nf > 1 cor- 
rections 

The differential Bhabha-scattering cross sec- 
tion with respect to the solid angle f2 can be writ- 
ten by means of an expansion in the fine-structure 
constant a, 
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where <jq is the Born contribution and er^ (i 
1.2 represent the higher-order radiative cor- 
rections. If we are interested in the NNLO virtual 
contributions, we need to evaluate the diagrams 
of Figure 2, where the fermion self-energy is re- 
quired for wave-function renormalization. Note 
that results for the photonic vacuum polarization 
diagrams can be found in [53] . 

After interfering the two-loop amplitude with 
the tree-level one, summing over the spins of the 
final state and averaging over those of the initial 
state, we get a large number of integrals. We 
use the DiaGen/IdSolver 34J implementation of 
the Laporta-Remiddi algorithm [TU] in order to 
reduce all the needed Feynman integrals to a lim- 
ited set of Mis. Apart from products of one-loop 
integrals, we get the eight Mis of Figure [I] as 
already pointed out in [13] . The corresponding 
Feynman integrals with n propagators Di are de- 
fined as follows: 
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where 7 is the Euler-Mascheroni constant and we 
introduced the shorthand notations 



{v}n = Vx,...,V ni 
V a b...c = V a + V h + . . . + V c . 



(3) 



In contrast to [23] , we do not consider Mis with 
scalar products in the numerators. We have then 
to allow for higher powers of propagators. Of 
course, there are algebraic relations between Mis 
with scalar products in the numerators and Mis 
with propagators raised to higher powers. 



We construct our MB representations using 
the standard approach described in [35]. The 
Feynman-parameter integrals are derived one af- 
ter the other for the two subloops. In each step we 
replace the sum over monomials in the Feynman 
parameters by an appropriate MB representation. 
Due to the relatively simple structure of the con- 
sidered diagrams, it is easier to begin with the 
propagator-type subloop. 

As far as the electron self-energy is concerned, 
we only need the sunrise topology with the elec- 
tron on its mass shell, pi = to 2 . A number of 
results for this mass configuration of the sunrise 
diagram can be found in the literature. Analytic 
expressions for the residues of the poles in di- 
mensional regularization and the finite parts were 
given already in [35) . The explicit result for the 
0(e) terms, where e = (4 — d)/2 and d are the 
space-time dimensions, can be found in |36j . Note 
that the inclusion of the 0(e) terms for the sun- 
rise Mis is mandatory when deriving the com- 
plete squared amplitude, since the reduction to 
Mis generates inverse powers of e. 

For the self energy, we use our MB represen- 
tation in order to reproduce the known result for 
the Mis. Its general form, for arbitrary powers Vi 
of the propagators, is given by 
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The integration contour is a straight line parallel 
to the imaginary axis separating the poles gener- 
ated by Ti and T4 from those coming from T 2 , 
T 3 , T 5 and r 6 . 



The two sunrise Mis are defined by the follow- 
ing values for the powers of the propagators, 

SE312Mlm = SE312Mlm(l, 1, 1), 
SE312Mlmd = SE312Mlm(l, 2, 1). (6) 

Having a MB representation at hand, one needs 
to perform an analytic continuation in e from a 
range where the integral is regular to the vicin- 
ity of the origin, uncovering the singular struc- 
ture on the way. This is done by an automatized 
procedure implemented in the Mathematica pack- 
age MB . m {37j . The resulting MB representations 
are verified numerically in the Euclidean region 
against the sector decomposition approach as de- 
scribed in [55] . 

For the sunrise Mis, a straightforward appli- 
cation of the Cauchy theorem to the MB repre- 
sentation of Eq. (J4j) leads to a sum over residua 
which can be easily evaluated. Therefore, we re- 
produced the results of In general, however, 
one has to deal with multiple MB representations. 
For the vertex and box Mis, the evaluation of the 
needed sums is far from being trivial. 

As explained in the introduction, our purpose is 
to calculate the integrals by assuming a hierarchy 
of all scales, namely m 2 <C M 2 <§; s, t. First of 
all we identify the leading contributions in the 
electron mass following the procedure described 
in [23]. Then, by using the Cauchy theorem to 
express the integrals through sums over residua, 
we evaluate these sums with the aid of XSUMMER 
[39]. 

For the sunrise Mis the results depend on one 
variable, R = m 2 /M 2 , and read as 
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For vertices, the external electrons are on their 
mass shell, p 2 — m 2 , i = 1,2, and we introduce 
the Mandelstam invariant s = (pi +P2) 2 (see Fig- 
ure [1} . Since each of the two vertices is related 
to two Mis, we have to consider 
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(9) 



We follow the same strategy employed for the 
sunrise diagrams. First of all we derive the exact 
multi-dimensional MB representation, and then 
we perform first a small-mass expansions in m s , 
and then in M 8 , defined as the ratios of the 
fermion masses and the centre-of-mass energy, 
m s = ~m 2 /s, M s = -M 2 /s. The MB repre- 
sentations read as 
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A careful analysis of the powers of m s and M s 
under the MB integrals leads to the following re- 
sults, 
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For box diagrams, the external momenta are 
again on their mass shell, and we have addition- 
ally (pi— P3) 2 = t. After introducing M t = —M 2 / 
t, the appropriate MB representation is given by 
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We have to compute two Mis, 

B512M2m = B512M2m(l, 1, 1, 1, 1), 

B512M2md = B512M2m(l, 2, 1, 1, 1), (20) 

and an expansion in the high-energy limit of the 
appropriate three-fold MB representations leads 
to the following results, 
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4. Summation techniques 

In any realistic computation we have to check 
the structure of the ultraviolet (UV) and infrared 
(IR) divergencies. By combining the Mcllin- 
Barnes method with recently developed summa- 
tion techniques we are able to evaluate exactly 
(i.e. without a high-energy approximation) the 
residues of the UV and IR poles for each MI. 

A simple example is enough to illustrate our 
procedure. We consider the following one-fold in- 
tegral in the complex plane, related to the single 
pole of V412Mlmd, 
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where we recall that M s = —M 2 /s, the integra- 
tion contour is a straight line parallel to the imag- 
inary axis, c = —1/2 and we introduced 
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After closing the integration contour to the right 
of the complex plane and taking residua, the in- 
tegral I can be written by means of two inverse 



binomial sums, 
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Inverse binomial sums were recently studied by 
means of the log-sine approach in [40]. Another 
approach was developed in [41] by generalizing 
the summation algorithms introduced in [42] . A 
straightforward application of these techniques 
leads to a compact result, 



I = 



1 - X M 



ln(y M )(l + 4M s )-ln(M s ), (26) 



1 + x M 

where we introduced the variables Xm and i/m, 

-s 



XM 



V4M 2 - s - 



V4M2 



V4M 



T2 



t 



VM 55 V4AP-t + V-f (27) 

As an example, after additionally introducing 
the following variables, 



%/4m 2 



V4m 2 - t - ^f~t 



V47 



t + 



(28) 



we get the non-approximated expressions for the 
residues of the poles of the two box diagrams de- 
fined in Eq. pi]) . 

^-ir-H(0;x m ), 



B-2 



m 2 1 



( 2 -2H(0,-l;x m ) 
2H(0;x m )\2H(-l;x, 

1 + 2/M 



1 - J/M 



H(0;y M ) 



(29) 



and 



5° 



x m Dm 



m 2 M 2 (1-0(1-1^) 



x ff(0;x m )JJ(0;2/ M ), 



(30) 



where we used the HPLs introduced in [T5] . 

For completeness, we add here also the exact 
expressions for the diveregent parts of the vertex 
Mis: 



1 f 1 + xm 



2 [ 1 — X M 
5^1 + x m 

2 1 ~ yy', 



H(0;x M ) + 2 + \nR 
■H(0;x m ). (31) 



5. Summary 

From [13] we know the table of Mis for mas- 
sive two-loop Bhabha scattering. We were able 
to express all the Feynman integrals occurring in 
the amplitude through these Mis by algebraic re- 
lations. We presented at the workshop all the 
planar two-loop box Mis. The Nf — 1 Mis has 
been published in the meantime [23] . In this con- 
tribution, we provide the expanded results for all 
the Mis entering the Bhabha-scattering ampli- 
tude with two fermion flavors in the limit of small 
fermion masses at fixed scattering angle. The Mis 
may be also found at our webpage [43] . These Mis 
were one of the last missing ingredients for the 
evaluation of the virtual two-loop contribution to 
the differential cross-section. 

The computation of the last nine non-planar 
two-loop box Mis is under way. 
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